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Scaling level-spacing distribution functions in the "bulk of the spectrum" in 
random matrix models N x N hermitian matrices and then going to the limit 
— > oo, leads to the Predholm determinant of the sine kernel sin7r(x — y)/7T{x — y). 
Similarly a scaling limit at the "edge of the spectrum" leads to the Airy kernel 
[M{x)k\'{y) - Ai'(x)Ai(y)] /{x - y). In this paper we derive analogues for this Airy 
kernel of the following properties of the sine kernel: the completely integrable system 
of P.D.E.'s found by Jimbo, Miwa, Mori and Sato; the expression, in the case of a 
single interval, of the Fredholm determinant in terms of a Painleve transcendent; the 
existence of a commuting differential operator; and the fact that this operator can 
be used in the derivation of asymptotics, for general n, of the probability that an 
interval contains precisely n eigenvalues. 
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I. INTRODUCTION AND SUMMARY OF RESULTS 



A. Introduction 

In this paper we present new results for the level spacing distribution functions obtained 
from scahng random matrix models of N x N hermitian matrices at the edge of the support 
of the (tree-level) eigenvalue densities when the parameters of the potential V are not "finely 
tuned." This universality class is already present in the Gaussian Unitary Ensemble. It is 
known that these distribution functions are expressible in terms of a Fredholm 

determinant of an integral operator K whose kernel involves Airy functions. 

It turns out that there are striking analogies between the properties of this Airy kernel 

^^^^^^ ^ A{x)A'{y)-A'{x)A{y) 
x-y 

where A{x) = \/XAi{x), and the sine kernel 

A sin7r(a; — y) 
71 X — y ' 

whose associated Fredholm determinant describes the classical level spacing distribution 
functions first studied by Wigner, Dyson, Mehta, and others [[7|,pO|,p^ . 



To be a little more explicit about how these kernels arise we remind the reader that in 
the GUE the probability density that the eigenvalues of a random N x N hermitian matrix 
lie in infinitesimal intervals about Xi, - ■ ■ ,Xn is given by 



where 



N-1 

KN{,x,y) = J2 Vk{,x)ipk{y) ■ 

k=0 



Here {ipk{x)} is the sequence obtained by orthonormalizing the sequence 

|a;''exp(-x^)| 



over (—00,00). More generally the n-point correlation function _R„(xi, ■ ■ ■ , which is 
the probability density that n of the eigenvalues (irrespective of order) lie in infinitesimal 
intervals about a;i, ■ ■ ■ , x„, is given by 

Rn{xi, ■■■,Xn) = det (K^ixi, a;j))ij^^... „ 

In particular Ri{x), the density of eigenvalues at x, equals Kn{x,x). 
Replacing K]^{x,y) by 

^ ( ^ ^ I ^ ^ 

-A AT I Z + ^ , , , Z + 



has the effect of making the point z in the spectrum the new origin and rescaling so that 
the eigenvalue density at this point becomes equal to 1. Now if z is fixed then 

Ri{z) = Kn{z,z) ~ -V2N 
vr 



as 00, and 



TT / Tlx ny \ 1 sinnlx — y) 

hm —^==Kn I z H 7==, z + ' — 



N^<^^2N \ V2N V2N J TT (x-y) 

so we obtain the sine kernel in this way. (This limit is obtained by applying the Christoffel- 
Darboux formula to Kn{x, y) and then using the known asymptotics oi ipN{x) as N ^ 00.) 
The "edge of the spectrum" corresponds to 2; ~ ±\^2N and one has there the scaling limit 
10,111 



1 / X /— y \ Ai(x)Ai(y)-Ai(x)Ai(y) 



We now present some of the analogies which we have found. (That such analogies 
between properties of the sine kernel and Airy kernel exist is perhaps not suprising since 
they are both scaling limits of the same family of kernels.) The first is the analogue of 



the completely integrable system of P.D.E.'s of Jimbo, Miwa, Mori, and Sato [|T7| when the 
underlying domain is a union of intervals. The second is the fact that in the case of the 
semi-infinite interval (s, 00) (the analogue of a single finite interval for the sine kernel) the 



Fredholm determinant is closely related to a Painleve transcendent of the second kind (the 
fifth transcendent arises for the sine kernel ^7^)- And the third is the existence of a second 
order differential operator commuting with the Airy operator K. (The existence of such a 
differential operator in the sine kernel case has been known for some time, see e.g. |T^, page 
201.) This last fact leads to an explicit asymptotic formula, as the interval (s, oo) expands, 
for the probability that it contains precisely n eigenvalues (n = 1, 2, . . .) (the analogue of 
results in p|,|3l||). Here are our results, stated in some detail. (A short announcement appears 
in m.) 



B. The System of Partial Differential Equations 

We set 

m 

= U ('^2i-i, a2j) 
i=i 

and write D{J; X) for the Fredholm determinant of K acting on J. We think of this as a 
function of a = (ai, . . . , a2m)- Then 

2m 

4 log D( J; A) = -J2i-iyR{aj,aj)daj (1.1) 

where R{x, y) is the kernel of the operator K{I~K)^^ and da denotes exterior differentiation 
with respect to the a^-'s. We introduce quantities 

g,- = (J - K)-^A{a,) , = (/ - K)-^A{a,) , (1.2) 



(which are the analogue of the quantities r±j of [|T^]; see also [|,|T^,^,^]) as well as two 
further quantities 



u 



(A, (/ - K)-'a) , = (a, (/ - K)-'A) (1.3) 
where the inner products refer to the domain J. Then the equations read 
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dqj_ 
dak 
9Pj 
dak 

daj 

daj 

du 
daj 
dv 
daj 



-1 



aj — ak 
,k (IjPk - PjQk 



E(-i 
E(-i 



aj — ak 
k QjPk - PjQk 



aj — ak 
,k QjPk - PjQk 



aj — ak 



Pk (j 7^ k), 

Qk + Pj - QjU , 
Pk + ajQj + pjU - 2qjV . 



-lyq' 



3 ' 



-lypjQj . 

Moreover the quantities R{aj,aj) appearing in are given by 

R{a„ a,) = ^(-1)'^ iWk-Pjqkf ^ ^2 _ ^^^2 _ 2^^.^^.^ ^ 2g|t; . 



k^j 



Ctj — Clk 



These equations are derived very much in the spirit of see also H,|l5|,^ 



:i.4) 
:i.5) 
:i.6) 

:i.7) 

:i.8) 
:i.9) 



(1.10) 



C. The Ordinary Differential Equations 

For the special case J = (s, oo) the above equations can be used to show that q{s; A) 
(the quantity q of the last section corresponding to the end-point s) satisfies 



q" = sq + 2q^ 



ds' 



:i.ii) 



with 



g(s; A) ~ -\/A Ai(s) as s ^ oo. 

This equation is a special case of the P// differential equation 
derive for R{s) := R{s, s), which in view of ( p..l|) equals 

d 



;i.l2) 



T^J23[| . One can similarly 



the third-order equation 



^logfl(J;A) 



2\r R' 



;i.i3) 



It is also the case that 

R'{s) = -q{s-X)' (1.14) 

and this gives the following simple formula for D{J; A) in terms of a P// transcendent: 

D{J; A) = exp - s)q{x; dx^ . (1.15) 

This is much simpler than the corresponding representation of D{J; A) for the sine kernel 
in terms of a Py transcendent. The fact that q{s; A) satisfies ( |1 . 1 1\) can also be obtained by 
combining some results in [|T],^,pr3|. Thus in this case of a semi-infinite interval our results 
connect with inverse scattering. 



D. Asymptotics and the Commuting Differential Operator 

Again we take J = (s, oo) and consider asymptotics as s — oo. (Asymptotics as 
s — > oo can be obtained trivially from the Neumann series for (/ — K)^^.) From the random 
matrix point of view the interesting quantities are 

(_l)n gn 

nl aA" A=i 

This is the probability that exactly n eigenvalues lie in J. Using our differential equations 



plus the known asymptotics of q{s; 1) as s — — oo [|T3| we obtain the asymptotic expansions 



^, 1 2 1 9 189 21663 , 

R(s;l) \ -H -tH (s^-oo). (1.16) 

^ ' ' 4 8 s 64 128 512 5^0 ^ ' ^ ' 

(where R{s; 1) denotes R{s) corresponding to A = 1) and 



E{0;s) = D{J;l)^j-^exp (s'/12 



^_s)l/8 

3 2025 2470825 
^' ~ 26^ + 213^" 219 s9 ^' 



■) is ^ -oo) (1.17) 



where tq is an undetermined constant. 

The analogue of this formula for the sine kernel was obtained by Dyson The ana- 
logue of our constant tq in the expansion was found to be equal to iC{s) is the 



Riemann zeta function). This constant was obtained by scaling a result of the second author 



on Toeplitz matrices. Such a result is not available to us now. Using ( |1.16| ) and the fact 
that R{s) is the logarithmic derivative of E{0] s), it is straightforward to derive an integral 
representation for logro. This integral can be evaluated approximately by numerically inte- 
grating ( |1.11D and using ( p..l4| ). This leads to the value logro — —0.136540. We conjecture 
that in fact 

ro = e^'("^) 2^ = 0.87237 14149 54127 .... 

This agrees with the numerically computed value, up to its degree of accuracy. 
For asymptotics of E{n; s) for general n we introduce 

E(n;s) 



r{n; s) :- 



E{0;s) 



Successive differentiation of ( |1.11|) with respect to A, plus the known asymptotics of q{s; 1), 



allows us to find asymptotic expansions for the quantities 



<ln[S) 



(1.18) 

A=l 



(9A" 

(for the analogue in the sine kernel case see 0); and these in turn can be used to find 
expansions for r(n;s). One drawback of this approach is that yet another undetermined 
constant factor enters the picture. (In [pj^ Toeplitz and Wiener-Hopf techniques, not 
available for the Airy kernel, fixed this constant.) Another drawback is that when one 
expresses the r{n] s) in terms of the g„(s) a large amount of cancellation takes place, with 
the result that even the first-order asymptotics of r{n; s) are out of reach by this method 
when n is large. 



There is, however, another approach (briefiy indicated in p| , and with details in ||2^ , for 
the sine kernel case). We have 

^.£. (l-A!:^^:^A.J 

where Aq > Ai > ■ ■ ■ are the eigenvalues of the integral operator K (with A = 1). Now just 
as the operator with the sine kernel commutes with the differential operator for the prolate 



spheroidal wave functions [|T4|, so does the Airy operator commutes with the differential 
operator L given by 



Lf{x) = ((x - s) fix))' - xix - s)f{x) . (1.20) 

An application of the WKB method, plus a trick allows us to derive the following asymptotic 
formula for Aj with i fixed: 



1 _ A, ~ ^ 25^+3 ^3^/2+3/4 (^-|t3/2 j ^ _2t ^ _oo). (1.21) 

(The analogue of this for the sine kernel is in |jll|.) From this and (|1.19|) we deduce 

l!2!...(n-U! . _ /8..,,,, 



-■(". -2*) ~ ;:;.2,^::,„,j exp [^nt"') , (1,22) 

(Note that this can be used to fix the constant mentioned in the last paragraph.) 



E. Probability Density for the n^^ Largest Eigenvalue 

An elementary probability argument shows that if F{n; s) {n = 0, 1, ■ ■ ■) denotes the 
probability density for the (scaled) position of the n*'^ largest eigenvalue (n = is the 
largest), then 

dE{0; s) 



F{0;s) 



ds 



and for n > 1 



^/ N N dE(n: s) 

F{n;s) - F{n-l,s) = ^ ' ' 



ds 

From the preceding expressions for E{n;s) in terms of q{s; 1) and qn{s) we can derive 
formulas for F(n; s) which can in turn be numerically evaluated to produce tables of values 
for these probability densities. We have done this for n = and n = 1 and in Fig. 1 we plot 
E{0; s) and E{1; s) and in Fig. 2 we plot F(0; s) and F(l; s). 
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II. THE SYSTEM OF PARTIAL DIFFERENTIAL EQUATIONS 



A. Derivation of the Equations 

A central role in the derivation of our equations is played by the commutator identity 

L, (J - K)-'] = (J - K)-' [L, K] (/ - K)-' (2.1) 

for any operators K and L. Taking L to be D (= d/dx) or M (=multiplication by the 
independent variable) give two of the three basic ingredients from which the equations are 
derived. The third is the equally simple fact that if K depends upon a parameter a then 

|-(,_,o-'^-(/-xr'^(/-Arv (2.2) 

It will be very useful for us to think of K as acting not on J but on (—00, 00) and to 
have kernel 

K{x,y)Xj{y) (2.3) 

where Xj is the characteristic function of J. We continue to denote the resolvent kernel of 
K by R{x, y) but note that while it is smooth in x it is discontinuous at y = aj. Thus the 
quantity R{aj,aj) appearing in ( p.. ID must be interpreted to mean 

lim R{aj,y); 

similarly for pj and qj in ( [I.2[ ) and like quantities appearing later. The definitions of u and 
V must be modified to read 



u 



{ax J, (/ - K)-' a), v= {AXj, (/ - K)-' A') . (2.4) 
Notice that since 

{I - K)-' A = {I - K)-' AXj in J 
this agrees with the original definitions of u and v given by ( |1.3|) . 



We can think of K as an operator taking smooth functions to smooth functions and 
so the operators KD and DK make sense, the former having distributional kernel. The 
transpose operator takes distributions to distributions. 

With these preliminaries out of the way we begin with the derivation of the representa- 
tions of the quantities i?(aj, a^) in terms of the g^, pj, u and v. We introduce the notation 



Q={I-Ky^A, P={I-Ky'A 



-1 M 



(2.5) 



so that qj = Q{aj), pj = P{aj). (For the moment we ignore the dependence of Q and P on 
the parameter a = (ai, ■ ■ ■ , ci2m); this will come later.) 
The first commutator relation we use is the simple 

[M,K] = {A{x)A\y)-A'{x)A{y))Xj{y) 



(where = means "has kernel equal to"). Hence from (Tl) with L = M 



M,{I-K)-' =Q{x)(l-K') A'xAy)-Pix)iI-K)-'AxAy) 



(2.6) 

(The transpose here arises from the general fact that if L = U{x)V{y), then T1LT2 = 
TiU{x) TgVd/).) If we use the fact that 



(^I - K') ^ AXj = {I - K)-^ A on J, 



(2.7) 



we deduce from (12.51) and (12.61) that UTS 



Q{x)Piy) - Pix)Qiy) ^ j , ^ 

R{x,y) = {x,y e J,x ^y) 



In particular, therefore. 



x-y 



R{aj, Ofc) = [J r k), 



(2.8) 



R{aj, aj) = Q'iaj)pj - P'(%)gj 



dx 



(2.9) 
(2.10) 



To compute Q'{x) and P'{x) we consider first the commutator [D,K]. If an operator L 
has distributional kernel L{x,y) then 
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It is an easy verification, using the equation A"{x) = xA{x), that 
and so (recall the definition ( p.3|) of the kernel of K) 



(2.11) 



(2.12) 



[D, K] = -A{x)A{y)Xjiy) - J2i-l)'K{x, a,)5{y - a,) . 



Hence, applying (|2.1| ) with L = D (and recalling definitions (|2.5| )), 



D, (/ - K)-^] = -Q{x) (/ - K') AXj{y) - ^(-l)^i?(x, a,)p(a,, y) (2.13) 



where 

p{x,y) = 5{x -y)+ R{x,y) 

is the distributional kernel of (/ — K)~^. 

We proceed with the computations of Q'{x) and P\x). Differentiating the first relation 
in (|2.5|) and using ( p.l3| ) give 



Q' = D{I - K)-^A = P+[d,{I- Xy'^] A 

= P-q((I- K') AX J, a) - Y,{-lfR{-, ak)qk 



Thus (recall (|J)) 



Q'{x) = P{x) - Q{x)u - Y.{-lfR{x, ak)qk 

k 



Similarly 



(2.14) 



P' = {I- KY^A" + \d, (/ - KY^] A' 



Since A"{x) = xA{x) this is equal to 



M (/ - KY^ A - \m, (J - KY^] A + \d, (/ - KY^ 



A' 
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We now use formula (pl6| ) and (|2.13|) to deduce that 



P'{x) = xQ{x) - 2Q{x)v + P{x)u - J2{-l)''R{x, ak)pk ■ 



(2.15) 



Setting X = aj in (|2.14| ) and (|2.15|) and using (|2.1CI|) give 



R{aj, ttj) =p]- a^q^- + 2q^-v - 2pjqjU + ^{-lfR{aj, ak)iqjPk - Pjqt) ■ 

k 

In view of ( p. 91) this is precisely equation ( |1.10|) . 

Most of the work is already done and we can derive the equations ( [I.4|) -( pr9D very quickly. 
First, we have the easy fact 



_d_ 
dai 



K 



-l)^K{x,ak)5{y - ak] 



and so by ( ^l2| 



dai 



{I-Ky' = {-lfR{x,ak)p{y-ak). 



(2.16) 



At this point we must keep in mind that Q and P are functions of a as well, and so we 
denote them now by Q{x,a) and P(x,a), respectively. We deduce immediately from ( |2.16| ) 



and the definitions (|2.5|) that 
d 



d 
— a) = {-l)''R{x, ak)qk, -^P{x, a) = {-lfR{x, ak)pk ■ 



dttk ' ' "v-j-^y-i^, 

Since qj = Q{aj, a) and Pj = P{aj, a) this gives 



(2.17) 



-l)''R{aj,ak)qk, 



dpi 



dak ^ ' Qa^^ 
In view of ( p. 91 ) these are equations ( |1.4| ) and ( |1.5| ). Moreover 



-\fR{aj,ak)pk U^k). 



dttj \dx dttj J ' ' 



x=aj ' dttj \dx dttj J 



and (|2A3), (^J^) and (pJfSD give 



dp 



— = ajqj - 2qjV + pju - '^ {-1)'' R{aj , ak)pk ■ 

j k^j 
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In view of ( p. 91) again, these are equations ( |1.6|) and (p..7|). 
Finally, using the definition of u in (p.4D , the fact 



and (|2.16| ) we find that 



^ = {~iyA{a,)q, + i-iy {ax J, R{; a,)) g, . 

But 

{AXj,R{-,aj)^ = J R{x, aj)A{x) dx = J R{aj, x)A{x) dx 
since R{x, y) = R{y, x) for x,y E J. Since R{y, x) = for x ^ J the last integral equals 

/oo 
R{aj, x)A{x) dx = A{aj) — qj . 
-oo 



Thus 

which is ( |1.8|) . Equation (|1.9|) is derived analogously. 

We end this section with two relations between u, v, the qj and the pj which in fact can 
be used to represent u and v directly in terms of the qj and pj. These relations are 



2v-u' = Y^{-iyql (2.18) 

j 

« = - E (P? - S^J - 2Pj?j^ + 2q^v) . (2.19) 
j 

To obtain the first of these observe that ( |1.4| ) and ( |1.6| ) imply 



(E^) (2.20) 



while from (lOD and (0 



^ {2v-u')=2i-iyq,{p,-q,u) . 
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Hence if we multiply both sides of (|2.20|) by 2{—iyqj and sum over j we obtain 

d 



It follows that the two sides of ( |2.18[ ) differ by a function of (oi, • ■ ■ , a2m) which is invariant 
under translation by any vector {s,---,s). Since, clearly, both sides tend to zero as all 
oo, their difference must be identically zero. 
To deduce (|2l9|) we add (|2l3|) and ( ^IBD to obtain 



(Observe that the kernel of [D, (/ — KY^] is {d/dx + d/dy)R{x, y).) Hence 



,a,) = -gJ = -(-iy 



du 
da A 



"3 

the last by ( |1.8| ). Multiplying by (—1)-' and summing over j gives 



From this we deduce, by an argument similar to the one at the end of the last paragraph, 
that 

E(-im«„«.) = -«. (2.21) 

j 

If we substitute formula ( |1.10|) into this we see that the resulting double sum vanishes and 



we are left with ( |2.19| ) 



B. Hamiltonian Structure 



The partial differential equations ( |1.4] )- (|1.9| ) can be rewritten in a compact form that 
reveals a symplectic structure with ( |1.10| ) defining a family of commuting Hamiltonians. To 
see this we first change the notation slightly 



I 



q2j = -2^2j, P2j = -ty2j, 
_ 1 

q2j-l - 2^2i-l, P2j-1 - y2j-i, 
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for j = 1, 2, . . . , m, and 



V = -Xo, u = yo, 



and introduce the canonical symplectic structure 

{xj,Xk} = {yj,yk} = 0, {xj,yk} = Sjk, 
for j, = 0, 1, . . . , m. Then equations (|L4|)-(|L9[) become 

daXj = {xj,uj{a)} and daVj = {yj, uj{a)} (j = 0, 1, . . . , m) (2.22) 



where 



and 



uj{a) = da log det (/ — K) 

2m 

= Y.G,{x,y)da, (2.23) 



Gj{x, y) = yj - -ajXj - Xjyjyo + -XjXo - t 2^ ' 



4 ^ ^ 4-^.- 4Z_. 

Furthermore, the G/s are in involution 

{G„Gk} = 0. 

This last result can be verified by a direct calculation, but as we discuss below, there is a 
better way to understand why such Hamiltonians are in involution. We can summarize the 
system of equations by saying, in words, that to find the Ofc-fiow of the coordinates xj and 
yj, one flows according to Hamilton's equation with Hamiltonian G^- The consistency of 
these equations (i.e. mixed partials are equal) follows immediately from this Hamiltonian 



formulation (see, e.g. p8[1). 



When equations ( |1.4| )-( pr9| ) are expressed in the Hamiltonian form ( |2.22| ), the equations 
are exactly analogous to the case of the sine kernel and in both cases the Hamiltonians 
are defined via ( p.23| ). For the JMMS equations (see, e.g. [^) this Hamiltonian system is 
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closely related to the Hamiltonian system studied by Moser |2J]. In fact in both cases the 
systems of partial differential equations can be formulated in terms of the classical i?-matrix 
approach to completely integrable systems. Such an approach identitfies the Lax operator 
and hence the commuting flow result above follows without any calculations. We refer the 
reader to |T2| for a discussion of this classical i?-matrix formulation of these systems of 
partial differential equations. 



III. THE ORDINARY DIFFERENTIAL EQUATIONS 

In this section we specialize to the case J = {s, oo) and derive the differential equations 
(|1.13|) and (|1.14|) . In the notation of the last section m = 1, Oi = s, 02 = 00. All 



terms in the formulas we derived which contain 02 vanish because of the decay of the Airy 
kernel at +00. Recall that we now denote the quantity q corresponding to the endpoint 
ai = s by q{s), and write R{s) for R{s, s). 

First, equation ( |1.14D is immediate from (|2.21|) , which now says that R = u, and ( p..8|) . 

For (|1.13|) we observe that (|L^) and (|1.7|) in the present case read 



q'=p-qu, (3.1) 
p' = sq + pu — 2qv , (3.2) 



whereas ( |1.10| ) reads 



Thus, by ( |LT1D , 



R = pq' — qp' , 



R' [qj • 

Taking the logarithmic derivative of both sides of ( p,.14| ) gives 

IR" _ q' 
2'RI ~ 7 

and by ( p.l|) this equals 

16 



(3.3) 



Thus 



P 



1 

2 



u 



V 



-R. 



n -R' 



By (U) this is just ( |LT^ ). 

Finally, to obtain ( |1.11| ) we differentiate and use ( 
find by an easy computation that 

q" = sq + + q{u^ - 2v) . 



3), (O) and (pp. We 



Conveniently, we have ( 2.18| ), which says in this case that — 2v = q^ . Thus ( |1.11 
established. 



IS 



IV. ASYMPTOTICS 
A. Asymptotics via Painleve 

Hastings and McLeod (see also PJ^) have shown that g(s; 1), the unique solution to 



equation ( |1 . 1 1| ) which is asymptotic to Ai(s) as s ^ +oo, is asymptotically equal to J—s/2 



as s — oo. In fact there is a complete asymptotic expansion in decreasing powers of —s, 
beginnning with this term |T9[ . Successive terms are easily computed. This result and later 



ones, become a little easier to state if we write s = —t/2: 

, ,^ ^ r(. 1 73 10657 13912277 1 ,\ , , 

Squaring and integrating and using (|1.14| ) we deduce 



, 1 2 1 9 189 21663 1 , , 

^^^) = 4^ +^-8l+647-T287+5T27^ + ^(;i^) ^ "^^ 



where c is a constant. Substituting this into ( |1.13| ) shows that c = 0. Thus we obtain (|1.16|) 
and then, after integrating and exponentiating, (|1.17|) . 
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To obtain the asymptotic expansions of Qnis) in ( |1.18|) , we differentiate ( |1 . 1 1| ) with 
respect to A and set A = 1. We obtain 



q'l - {s + 6go) gi = 



(where qo{s) = q{s; 1)). Using ( ^?T| ) we find two hnearly independent solutions of this equa- 
tion, one exponentially large at — oo and one exponentially small. Assuming the exponen- 
tially large one actually appears with a nonzero factor we obtain the asymptotic expansion 

exp(|t3/2) . ^7 ^5^3 850193 407117521 \ 
gi(-t/2) - ci -^^ [1 + + + 21034^9/2 + 21535^6 + ■ ■ 7 ^^-2^ 

where ci is a constant factor to be determined. 

Successive differentiations of ( |1.11|) with respect to A, followed by setting A = 1, yields a 



sequence of linear differential equations for the 

where is a polynomial in the qm with m < n. For n > 1 a, particular solution dominates 
all solutions of the homogeneous equation. Therefore, in the asymptotics, no new unknown 
constant factors are introduced (unless we want to go "beyond all orders"). 

To deduce expansions for the r{n; s) we must differentiate ( |1.15| ) n times with respect to 
A and set A = 1. We find expressions of the form 

r{n;s)= {x - s)i){qo,- ■ ■ ,qn) dx (4.3) 

J s 

where ipn is a polynomial in go; ■ " " ; In- In particular we obtain (after setting s = —t/2) 

r(l; -t/2) = ^ r (t - x)qo{-x/2)qi{-x/2) dx . 



2 



00 



Substituting ( |4.1| ) and ( |4.2| ) into this integral gives the asymptotic expansion for r(l; —t/2). 
The first order result is 

exp(t3/2/3) 
r(l;-t/2) ~Ci-^^^5^ 



and comparing this with (|1.22| ) reveals that 

1 

Cl 



2V27r 

We now state the expansions we have obtained in this way with the help of Mathematica. 
For Qn with > 2, as t ^ +00 

n\ exp(f ^3/2) / 59^ _ 36 

gn(-V^) - 23n/2^n/2 ^3n/4-l/2 1''-^ ^ ^ + 

3481n2 + 8244n - 8496 , n/2^^ .^ 
__ t-3 + or^/^) j . (4.4) 

For n = 1, 2 the leading order term is as above. For n = 1 the correction terms are given 
above in ( [4.2|) . For n = 2 the coefficient of is as given above but the coefficient of 
is 5461/288. (Similar exceptional values of n occurred in analogous expansions for the sine 
kernel p[.) These expansions have been computed for n < 10 and are undoubtedly correct 
for all n. (In the analogous expansions were proved for all by a backwards induction 
argument — such arguments should extend to the present case of g„.) 

For r{n] s) we have obtained the following sharpening of ( |1.22|) : as t — > +oo 
_ l!2!---(n-l)! exp(tt3/2) / n(34n^ + 31) /, 

n^(1156n^ + 6608n^ + 11509) ^ q^^-9/2)^ _ (4 5) 

This was established for n <Q. As was mentioned in the Introduction, when we substitute 
expansions ( [4. 4] ) into ( |4.3| ) a large amount of cancellation takes place, so that r{n] s) is much 
smaller than might be anticipated. This also means that to deduce even the ffist-order 
asymptotics in ( [4.5|) for moderate values of n one has to go very far out in the expansions 
( [4.4|) . Thus it is important that we have an alternative approach which yields asymptotic 
results for all values of n. 



B. Asymptotics via the Commuting Differential Operator 

We begin with the integral identity 
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K{x,y)= A{x + z)A{y + z)dz. (4.6) 
Jo 

This appears in but we give a simple proof here. Think of ( |2.12| ) as a differential equation 
for the "unknown" function K{x,y). Since 



one solution is the right side of (|4.6|) and, of course, another solution is the Airy kernel 
K{x,y). The difference is therefore of the form (p{x — y) for some function (p. Since both 
sides of ( ^l6|) tend to zero as x and y tend to +oo independently, we must have (p = 0. 
Rewriting (^l6| ) as 



/■oo 

K[x,y) = / 74(x + z — + 2/ — s) (i^; 

J s 



shows that K as an operator on (s, oo) is the square of the operator with kernel A{x + y — s). 

Integration by parts shows that an integral operator with kernel L{x,y) on (s, oo) com- 
mutes with a differential operator 

^«(a:)A + ^(^) 

if a{s) =0 and if 

d'^L dL d'^L dL 

"^^^^W ^ "^'^^^dy ^ ^^^^^ " "^"^^ + "'(^) + ^^^)^ • 

(Of course we also require, in the end, appropriate vanishing at oo.) If we set L{x,y) = 
A{x + y — s) and use the fact A"{x) = xA{x) we can easily check that a solution is given by 

a{x) = X — s, (3{x) = —x{x — s) 

and so we obtain the commuting differential operator ( |1.20| ). 



Here is an outline of how ( |1.21| ) is found. If we let /i, /2, ■ ■ ■ denote the eigenfunctions 
of L then the logarithmic derivative of each Aj (thought of as a function of s) is expressible 
very simply in terms of the corresponding /j. The asymptotics of these eigenfunctions, and 
so of the logarithmic derivative, is established by WKB techniques. Integrating, and using 

20 



the fact that each Aj 1 as s — — oo, gives (|1.21| ). The WKB argument we give is quite 
heuristic although it could very likely be made rigorous. 

We make a preliminary change of variables, replacing x hj x + s (so that L acts on 
(0, oo)) and s by —2t. The eigenvalue problem can then be written as 



{xf'{x))'+{f,-{x-ty)fix) = 



(4.7) 



or in the equivalent form 

xg"{x) + 



/i 



(4.8) 



The eigenvalues satisfy < /io < < ■ ■ ■ and oscillation considerations applied to ( |4.8| ) 
show that /ij = 0(t^/^), for each i as t oo. So we think of /i in our equations as 0(t^/^), 
we normalize our eigenfunctions so that /(O) = 1, and proceed to find asymptotics for large 
t. 



1. The region x « t 

If we make the substitutions 

fix) = h{t'x) = h{y) 

then ( |4.7| ) becomes 

{yh'{y))'-h{y)=t-' (^^ - 2| - /.^ %) . 

The solution of the equation with the right hand side replaced by which satisfies h{0) = 1 
is /o(2y/y) where Jq is the usual modified Bessel function. The actual solution h{y) will be 
asymptotic to this in any interval over which the integral of the factor of h{y) on the right 
hand side of the equation is o(l). Recalling that /i = 0(t^/^) we see that y « t^^"^ suffices. 
Hence 

f{x) ~ Jo (2tv^) , (x « t-i/2) (4.9) 



21 



2. The region x » t , t — x » t 



1/4 



If in equation ( |4.8| ) we set 



a(x) = ^ —, y = f \fa{z) dz , hiy) = aixV^'^gix) 

X Jo ^ 



then the equation for h{y) is 

h"{y) - h{y) 



la" 5 a'^ 1 
4 16 Ax'^a 



h{y) 



where on the right side ' = Since dy = q;(x)^/^ dx, the integral with respect to y (of the 
factor of h{y) on the right side of the equation) over any interval equals 



1 a" 5 a 



l2 



1 



dx 



over the corresponding x interval. This will be o(l) as long as, in this interval, a; >> 
and t - X >> t^/'^. 

In this same region y » 1 and so 

h{y) ~ ae^ (4.10) 

for some constant a = a{t). In case x << 1 we have y = 2t^/x + o(l) and so 

f(x) = x-^'^g{x) ~ at-^'^x-^'^e"'^ , (r^ « x « 1). 

Comparing this result with (|4.9| ) in the overlapping region << x << t~^/^ and using the 
known asymptotics of Io{z) as z ^ oo shows that a ~ 1/2-^/??. Hence from 10|) 

1 (t-x)-i/2 



(x >> t~^,t-x » t^/^) 



We also used here once again the estimate /i = 0(t^/^). Because of this same estimate we 
can write the integral in the exponential as 



[ 2{t-z)^ \{t-z) 
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Because t — x >> t^^^ the last term in brackets contributes o(l) to the integral, and so we 
obtain 



3 2 Jo ^{t - z) 

Replacing z by tz^ shows that the integral here equals 

log V 

t — X 

Thus 

m ~ ^ '[Zn, - \^"'] (4-11) 

{x » r'^,t-x» t^/^). 

Let us see what happens if x is not far from t, if 

t'" «t-x«t"\ 

Then since 

2teV2 _ 2^3/2 ^ 4 4/3 _ _ ^ 

3 3 2^ ^ ^ ^ 



in this region, we find 



5. r/ie region x — t » 

^1/4 



This is quite similar to the preceding. The main difference is that instead of making the 
variable change y = J a^xY^'^ dx in ( |4.^ ) we now set 

2 



y = f a;3/2 _ 2^2,1/2 _ r 



a U - z''^ + tz 



1/2 I ^^-1/2 
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Since an eigenfunction cannot grow exponentially at oo only the negative exponential e ^ 
can appear in the asymptotics and we find that 



m ~ 6^^^^exp htx'/' - Ix'/' + 



for some constant b = b(t). Now the integral in the exponential is 

. 2 



+ Vi 



X — t 



+ o(l) 



and so 



fix) ~ b- 



ix-t)^'-'^'''-^ 



exp 



f{x) ^b2-''""H'''/^ 



X — t 



exp |-r'" - -(x - t)' 



dz 



and specializing to the region t^/^ << x — t « t^l'^ gives 

{t'/^«x-t«f/') 



(4.13) 



(4.14) 



4- The region \x — t\ « t 



If in equation (|47q ) we set 



g{x) = h {V2t-^'\t - x)) = h{y) 



and write e = {\/2t^^^) ^ then the equation becomes 



(l-ey) h"{y) + 



y 



h{y)=Q. 



2ti/2 4 A{l-ey) 

Now we expect that for bounded y, and therefore also if \y\ oo at some rate (depending 
upon e), the solution of this will be asymptotic to a solution of 



Weber's equation. Now unless 
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^ -^ + 1 (^ = o,i,---) 



2ti/2 2 

there is no solution which vanishes at both ±00. But ( [4. 121 ) and ( |4.14| ) show that in their 
regions of vahdity (which overlap with our present region) f{x) is a factor depending on 
t times a factor tending exponentially toOasx — t— s>±oo. Hence our eigenvalues must 
satisfy 

-^ + l + o{l) 



2tV2 2 

and h{y) is asymptotically a constant times the parabolic cylinder function Di{y) (see, e.g. 
0, Chp. 8). Thus (since x ~ t in the present situation) 

fix) ^ cDi[V2r^/\t-x)) (4.15) 

for some constant c = c{t). Using the known asymptotics of Di we find 

fix) ~ c [V2t-^/\t - x))'exp |-ir^/2(t - x)2| (4.16) 



Comparing this with ( 4.12| ) gives 

2-1*-^ 3. 3 /4 



c 



rt*-iexp -t3/2 . (4.17) 



We also find, using (|4l^), that 6 ~ 1/2 v^. 



5. The asymptotics of Aj 



If we look at the asymptotics of fi{x) (the eigenfunction associated with the eigenvalue 
yUj of L, normalized so that /i(0) = 1) given by (^l9|), ( |4.11| ), ( [4.13|) , and ( [4.16|) we see that 
the main contribution to /q°° fi{x)'^ dx comes from any region 



t-^l^\t-x\ < r]{t) 



as long as ■r](t) ^ 00. For some such region we have, by ( [4. 15] ) and ( |4.17| ), 



f^ix) 



2-i 



i-2 



X 
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Using this, and the fact Jf^ Di{x) dx = vSvri!, we deduce 

r U{xfdx ~ 4=2-'^-^^^-^ exp f-t^/^) . (4.18) 

^0 VTT \3 / 



This was the goal of all that went before. 

Since L has simple eigenvalues and commutes with our integral operator K (translated 
to (0, oo)) the set {/j} is the set of eigenfunctions of K corresponding to its eigenvalues {Aj}. 
However /j corresponds to /Xj and there is no assurance that the corresponding eigenvalue 
of K is its fi^ largest, which we have denoted by Aj. We shall eventually show that this is 
actually the case. The first lemma of this section shows that the eigenvalues of K are simple, 
and this is crucial for the argument we shall later use. (The proof of the corresponding result 
for the sine kernel is in 



Lemma 1 // /i and f2 are eigenfunctions of L corresponding to distinct eigenvalues, then 
they are eigenfunctions of K corresponding to distinct eigenvalues. 

Proof. The eigenvalues of K are the squares of the eigenvalues of the operator on (0, oo) 
with kernel A{x + y — 2t). So we have to show that if for some i' we have either 

A{x + y-2t)f,{y)dy = uf,{x) {t = 1,2) (4.19) 



or 

/•oo /"OO 

/ Aix + y-2t)f^iy) = lyf^ix), / Aix + y - 2t)f2iy) dy = -uf2ix) (4-20) 

Jo Jo 

then fi = f2- (Recall that all eigenfunctions are normalized to satisfy /(O) = 1.) 

Assuming first that ( [4.19| ) holds, apply (ff' jdx^ to both sides of the identity with i = \ 



and integrate the resulting integral by parts twice. What results is 

("OO 

vf'l{x) = -A\x - 2t) + A{x - 2t)f'M + A{x + y- 2t)f'^{y) dy . 



If we multiply both sides by f2{x) and integrate we obtain, using (|4.19| ) with i = 2 (and its 
differentiated version), 
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Jo Jo 
Thus, since u ^ (this follows easily from the fact that the Fourier transform of A{x) is 
nonzero), we have /((O) = f^i^)- But it is clear from (|4.7] ) that for any eigenfunction / we 
have /'(O) =t^ — ^. Thus /i and /2 correspond to the same eigenvalue of L and so are equal. 
The case when ( 4.2CI| ) holds is even easier. Differentiating both sides of the first relation 



once and integrating by parts gives 

roo 

vf[{x) = -A{x-2t)- / A{x + y-2t)f[{y)dy. 

Jo 



Multiplying both sides of this by /2(a;) and integrating, using the second relation of ( |4.20D , 
we obtain 



oo POO 

f'i{x)f2{x)dx = -u + u f2{y)f[{y)dy , 
Jo 



contradicting i/ 7^ 0. ■ 

Since, as we now know, the eigenvalues of K are simple, there is a permutation a of 
N = {0, 1, 2, ■ ■ ■} such that /j is the eigenfunction of K corresponding to its a{iy^ largest 
eigenvalue \a{i)] this permutation is independent of t or else, by the continuity of the eigen- 
values in t, there would be a multiple eigenvalue for some t. 

Up to now our kernel K corresponded to an arbitrary factor A. Of course formula (|1.21|) 
refers to the case A = 1. 

Lemma 2 For each i we have Aj — * 1 as s — — 00. 
Proof. In this case A = 1 formula f|4.6|) reads 



roc 

K{x,y)= Ai{x + z)Ai{z + y)dz. (4.21) 
Jo 



For the purposes of this proof we denote by A (resp. K) the operator on (— 00, cxd) with 
kernel Ai(x + y) (resp. K{x,y)) and by Ps the projection from L2{—oo, 00) to ^2(5, 00). So 
the Aj are the eigenvalues of PsKPg. Now it is known [|T3| that A^ = I, and ( ^.211 ) says 



that K = APqA. Thus K"^ = APqA^PqA = APqA = K. In other words, is a projection 
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operator. (This is not surprising since its kernel is a scaling limit of the kernels KN{x,y) of 
projection operators.) Clearly it has infinite-dimensional range. It now follows easily from 
the minimax characterization of the eigenvalues that, for each i, the i^^ largest eigenvalue of 
PgKPs tends to 1 as s ^ — oo. ■ 

We now state a lemma which is the basis of a trick used by Fuchs [11| in his derivation 
of the asymptotics of the eigenvalues of the sine kernel. 

Lemma 3 Let J{x, y) he a symmetric kernel, on a fixed interval, depending smoothly on a 
parameter t. Let g be an eigenfunction normalized so that J g{xY dx = 1 and let A be the 
corresponding eigenvalue. Then, with the subscript t denoting d / dt , we have 

At = / / Jt{x,y)g{y)g{x)dydx. 



Proof. From 

J J{x,y)g{y)dy = Xg{x) (4.22) 

we get 



Multiplying both sides by g{x) and integrating over x, we obtain (using the normalization 
of g and ( ^1221) ) 



Jt{,x,y)g{y)g{x) dydx + / / J{x,y)gt{y)g{x) dxdy = Xt + X J gt{x)g{x) dx 

= Xt+ / J{x,y)g{y)gt{x)dydx. 



By symmetry of J the two double integrals involving it cancel, and we obtain the statement 
of the lemma. ■ 

We apply the lemma to the opertor with kernel A{x + y — 2t) on (0, oo) whose square 
is our Airy operator (translated to act on (0, oo)). Thus we may write its eigenvalues as 



Ao-(j), where the square roots might have either sign. The lemma gives 
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oo roo 



JQ 



A'{x + y- 2t)gi{y)gi{x) dydx 



(4.23) 



where 



coo 'I 1/2 

g,{x) = U{x)/\ 1^ h{yfdy\ 



But from 



A{x + y- 2t)gi{y) dy = \[K(i)gi{x) 
we obtain upon differentiation that 

A'{x + y- 2t)gi{y) dy = ^)^)g[{x) 

so ( ^.231 ) may be rewritten 



\l K{i) = -2^\„(^i) j g'i{x)g,{x) dx 



d_ 

The integral is of course —^gi{Oy. Recalhng that /j(0) = 1 we see that we have estabhshed 
the relation 

d 



dx 



or 



- log A.(,) = 2/11^" /,(a;fdx] 

We now use ([4.18| ) and integrate the resulting relation with respect to t from t to cx). 
Recalling Lemma ^ we deduce 

-logA.(.)~^2^^+='ti^+iexp f-V^^ 



II - \ 3 

precisely the right side of ( |1.21D . 

All that remains now is to show that cr(z) = i for all i. But it is clear from this asymptotic 
relation that i < j implies Ao-(i) > Xa{j) for large t (and so for all t), so that < cr{j). 
Thus a is order preserving, and since cr : N ^ N is onto we must have a{i) = i for all i. 
This concludes the derivation of ( |1.21[ ) . The deduction of ( |1.22| ) from ( |1.21D is completely 



analogous to the corresponding deduction for the sine kernel in ref. p8[, Sec. VIIB. 
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FIGURES 



FIG. 1. The probabilities E{0;s) and E{l;s). Of course, E{0;s) — 1 as s — +00. 

FIG. 2. The probabihty densities F{0;s) and F{l;s). The expected position of the largest 
eigenvalue is approximately —1.7711 with variance 0.8132. The expected position of the next-largest 
eigenvalue is approximately —3.6754 with variance 0.5405. 
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